On classification of discrete, scalar-valued Poisson brackets 



E. PARODI 



Abstract. We address the problem of classifying discrete differential- geometric 
Poisson brackets (dDGPBs) of any fixed order on target space of dimension 1. It 
is proved that these Poisson brackets (PBs) are in one-to-one correspondence with 
the intersection points of certain projective hypersurfaces. In addition, they can 
be reduced to cubic PB of standard Volterra lattice by discrete Miura-type trans- 
formations. Finally, improving a consolidation lattice procedure, we obtain new 
families of non-degenerate, vector-valued and first order dDGPBs, which can be 
considered in the framework of admissible Lie-Poisson group theory. 

Keywords: Discrete Poisson brackets, discrete Miura transformations, Lie-Poisson 
groups. 
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1. Introduction 

In this paper we deal with the following class of local Poisson brackets (PBs) 

{<)«n+fc}M = 9k{^n,---,Un+k), 0<k<M, 

defined on the phase space of infinite sequences 

u: Z — V 

with values in the target manifold of dimension N . The integer number M, 
called the order of the PB, can be seen as the locality radius, i.e. the radius of 
the maximum local interaction between neighboring lattice variables. These PBs 
have been introduced by B. Dubrovin in [3] (see also A. Ya. Mal'tsev, [10]), as a 
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discretization of the differential geometric Poisson brackets (DGPBs), defined on the 
loop space C{M'^) = M'^} by the formula 



M 



1.2) {u^ix),u^iy)}^^ = 5^ <7^^(u(x),u.(x), . . . ,u('\x))6^''^-'\x - y), 



k=0 



where i, j = 1, . . . , N and the functions g^j^ are graded-homogeneous polynomials. See 
for details the papers by S. Novikov and B. Dubrovin [4], [5]. 

PBs of type (1.1) are associated with lattice Hamiltonian equations of the following 
form 



m&L p=l "^"^ 

where if[u] = ^^g^ ^(^mi • • ■ ) u^+k) for some integer X > and the function h is 
defined on a finite interval of the lattice. In addition, we define the formal variational 
derivative as 

^ = ^ (1 + + . . . + T"^) Mun, . . . , u.+x) , 
where T is the standard shift operator, satisfying 

(1.3) T''/l(u„, . . . , Wu+k) = h{Un+r, ■ ■ ■ 

for any integer r. The local Poisson structures of many fundamental integrable sys- 
tems, such as the Volterra lattices, the Toda lattices, the Bogoyavlensky lattices (see 
Yu.B. Suris, [12]) belong to the class (1.1). However, the theory of such discrete PBs 
is much less developed than the corresponding of DGPBs (1.2) (see the survey of O. 
Mokhov [11] and references therein). 

A classification of first order (M = 1) PBs (1.1) has been provided in [3], whereas it 
seems that the higher order PBs have not been studied yet. Moreover, to the best of 
our knowledge, in the literature there are no examples of PBs (1.1) of order M > 2. 
In [3], a correspondence between the following first order PBs 

Q4^ {<><+l}i= C/l^(Un,U„+i) 

and certain Lie-Poisson groups was discovered. More precisely, if the matrix g^^ is 
non-singular (i.e. det (^^■'(u^, u„_|_i) ^ 0) the PBs (1.4) are induced by admissible Lie- 
Poisson group structures on the target manifold Ai'^ (see Definition 5.3 below). 
Performing a consolidation lattice procedure, which is obtained by defining new varia- 
bles of a larger target manifold by the formulas v^^^ = ul^M+p^ p = 0, . . . , M — 1, one 
can reduce any PB (1.1) to the form (1.4). However, this procedure leaves some 
unsolved questions: 

(i) what are the relations between PBs (1.1) of order M > 1 and admissible Lie- 
Poisson groups associated to their consolidations? 



(a) how to produce examples of such admissible Lie-Poisson groups? 

In the present paper, in order to give some partial answers, we classify scalar-valued 
(N = 1) PBs (1.1) of any positive order M, 

(1.5) {Un, Un+k}M = 9k{Un, • • • , Un+k), 1 < k < M. 

First, we observe that PBs of type (1.1) are invariant under local change of variable 

(1.6) u\^vi = v\un), z = l,...,iV, 
where the coefficients g^/^ transform according to the formula 

5'fe (Un, . . . , Un+k) I > ^2 (""■)5'r(^n, • • • , U„+fc)— -g (u„+fc). 

p,q=l ^""^ ^'^n+k 

Two local PBs will be therefore considered equivalent if they can be related by a 
change of coordinates of type (1.6). 

Let us be more precise about the classification of scalar- valued PBs (1.5). The co- 
efficients gk{un, ■ ■ ■ ,Un+k) satisfy the system of bi-linear PDEs imposed by the 
Jacobi identity (see Section 2). It turns out that any PB (1.5) is characterized by his 
leading order function gAiiun, ■ ■ ■ lUn+u)-, according to the following results. 

Lemma 1.1. For any PB of the form (1.5), there exist a set of coordinates (canonical 
coordinates) and an integer a > 0, such that the leading order reduces to the form 

gniUn, • • • , Un+Al) = f^{Un+a, • • • , «„+«+?), ^ = M - 2a, 
where the function f^ is either constant or given by the formula 

fn = fHun,---, Mn+c) = exp (Zn) . 

Here, 

(1-7) Zn = Tj Un+i , 

i=0 

and the parameters Ti, z = 0, . . . , satisfy a system of homogeneous polynomial equa- 
tions (see below Theorem 2.7). 

In the case when the leading coefficient gM{un, ■ ■ ■ , Un+Ni) is constant, it is not difficult 
to prove that all the other coefficients {5'fc}fc=i,...,M-i are also constant. The general, 
non-constant case is described by the following 

Theorem 1.2. The coefficients g^ of a non-constant PB (1.5) are given, in the 
canonical coordinates, by suitable linear combinations of the shifted generating func- 
tion f^, according to the following formula 

(m\n{a+p,a) \ 
KT'\ f^{Un,...,Un+d: 
s=max(0,p) / 
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where p = —a, . . . ,a, and the constants A 's can be expressed explicitly in terms of the 
parameters r's (see equation (2.8) , below). 

This Theorem provides a complete classification of PBs of type (1.5). Note that 
the functional form of the coefficients is fixed by the choice of a finite number of 
parameters r. 

In addition to the above results, we prove a Darboux-type theorem for PBs (1.5). 
This is done by considering the change of variable (1.7), that is a generalization of 
the local one (1.6). Notice that the formula (1.7) can be thought as a discrete analogue 
of Miura transformations, studied in the Hamiltonian PDEs theory (see [6]). 
Splitting all the variables into a + C, families according to 

(1-9) vl^^ = Z(a+0{n-i)+p, p=l,...,a + ^, 

by direct computation, we obtain that any non-constant PB (1.5) in the ^-coordinates 
can be reduced to the following simple form 

(p) (p) 1 / (p) \ 

,<+2r = ^OT-^exp^+i) 



„,(p) . _ _ _ 



exp(wi^^) + exp(t4+i) 



that are a + ^ copies of cubic Volterra PB (see Theorem 3.2 below). 
We consider next compatible pairs of PBs of type (1.5). Recall that a pair of PBs 
(Pi, P2) is said to be compatible (or to form a pencil of PBs) if any linear combination 
with constant coefficients ^Pi + VP2 is also a PB. This notion, first mentioned by 
F. Magri [9] and extended by I. Gel'fand and I. Dorfman [2] (see also [1]) provides a 
fundamental device for the integrability of Hamiltonian equations. In our setting, the 
study of compatible pair of PBs (1.5) might lead to the classification of bi-Hamiltonian 
lattice equations of type 

iln = F{Un-S, ■ ■ ■ ,Un, ■ ■ ■ ,Un+s) = {Un, Hi[u]} j^^^ = {Un, H2[u\} jy,^^ 

for some integer S > and local PBs of order Mi, M2. It seems that higher order 
lattice equations have not been studied yet except for the Volterra type equations 
(i.e. S = 1), analyzed by R.I. Yamilov and collaborators using the master symme- 
tries approach (see review article [15] and references therein). The following result 
describes some necessary conditions for the classification of pencil of PBs (1.5). We 
expect these conditions to be also sufficient. 

Theorem 1.3. Let us consider a pair (P, P') of non-constant PBs of type (1.5) and 
order M = 2a + ^ and M' = 2a' + ^' respectively, with M > M' . Then P and P' 
form a pencil of PBs only if a = a' and there exist coordinates on the manifold 
such that 



f^{Un+a, • • • , Un+a+i) = C^A^exp (^Yli=0 
f'^ {Un+a, . . . , Un+a+e) = 0"A/'exp (^XlLo '^i 



gM[Un, • • • , Un+M) = [Un+a, • • • , = ^A^^Xp ( ^^J^p Tj U^+i 

dM'iUn, . . . , Un+M') 
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where aM, o,"^^ some constants and Tp = = r^^^/^p, for any p = 0, . 

The structure of the present paper is the following: in Section 2 we describe our 
classification procedure, proving formula (1.8) for the coefficients gk- In Section 3, 
using the change of coordinates (1.7), suggested by the leading coefficient g^j, we prove 
that any PB of type (1.1) can be reduced to the cubic PB of Volterra lattice (Theorem 
3.2). In Section 4, we present some necessary conditions for the compatibility of 
(a, (^) -brackets (Theorem 4.2 and Lemma 4.3). Finally, with Section 5, we introduce 
the concept of admissible Lie-Poisson groups and describe how to produce some new 
classes of non-degenerate vector- valued (A^ > 1) first order dDGPBs (Theorem 5.6). 



2. Classification of {a, (^)-brackets 

This Section is devoted to the classification of the scalar-valued PBs given by the 
formulas 

(2.1) {Un,Un+k}M = 9k{Un, ■ ■ ■ , Un+k) , l<k<M. 

The coefficients gk{un., ■ ■ ■ ,Un+k) are locally analytic functions (see R. Yamilov [15] 
for the precise definition), satisfying, for all values of the independent variables Un, 
n e Z, the following bi-linear PDEs given by Jacobi identity 

[p, q] {{Un,Un+p} Un+p+q}} + {{ '^m ^n+p+g} i ^n+p} ) 

explicitly, 

X/i=0 dpi'^rn • • • ) Un+p) ,Un+p-i9q+ii'^n+p-i-i • • • j Un+p+q)''f 
~ X]i=0 dqiV'n+p-i • • • ) Un+p+q) ,Un+p+i9p+i{!^n-, ■ ■ ■ i Un+p+i) 

II 

'^^1=0 9p+qiUm • • • ) Un+p+q) ^u„+i9p—i{Un+ii ■ ■ ■ -j Un+p)-\~ 
~ Y2i=0 9p+q{U'ni • • • j Un+p+q) ,Un+p+i9i{U'n+pi ■ ■ ■ i Un+p+i) 

where p,q = 1, . . . , M and go{-) = 0, gk{-) = 0, k > M. 
Example 2.1. [Volterra lattice or discrete KdV equation] 

The well-known Volterra lattice (VL) [14] is defined by the following equations of 
motion 

iln = Un (Un+l - ti„-l) n G Z. 

Performing the local change of variable i — y log Un, we obtain equations 

(2.2) iln = exp(M„+i) - exp(M„_i), 

that admit the following bi-hamiltonian representation (see [7]) 

iln = {Un,Hi[u]}^ = {Un,H2[u]}^, 
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where 



^l(^i) = Efcexp(Mfc) {Un,Un+l}i = 1, 



(2.3) i„ {Un,Un+2}2 = exp(M„+i) 

^ {M„,u„+i|2 = exp(M„) + exp(M„+i). 

These PBs belong to the class (2.1). Moreover, the coefficients of the cubic PB { ■ , - ^^ 
are characterized by suitable linear combination of the leading order term, given by 
the exponential functions = /"(wn) = exp(u„). This will be the typical behavior 
of non-constant PBs (2.1). 

2.1. The leading-order coefficient. By the definition of locality radius, the leading- 
order term gu might depend on variables u„, . . . , Un+M, i-e- 9m = QAiiun, ■ • ■ i "^^n+Af)- 

Lemma 2.2. There exist canonical coordinates and an integer a > such that the 
leading term gM reduce to the form 

gAliUn, . . . , Un+Al) = fn+a = f{^n+a, Un+a+i) , ^ = M - 2a, 

where 

(2.4) M ^ ^ 

if f^ is a non-constant function. 

Proof. Let us consider the bi-linear PDEs [p, q], with p,q = 1, . . . , M. At ffist, we 
focus our attention on equation 

[M, M] {{Un, Uu+m} , Un+2M} = {Un, {Un+M, Un+2M}} 

that provides us 

log5'Af('U„, . . . ,Un+M),u„+M — log 9'A/(lin+Mi ■ ■ ■ ,Un+2M) ,u„+m — ^^(^n+Af) 

for some arbitrary function a{un+M)- Solving this logarithmic equation, we obtain 
the following factorization 

gAliUn, Un+Aj) = a{Un) f^^~"~'^ (Un+a, • • • , Mn+M-/3)a(M„+M) 

where loga(ii„),„„ = a(n„) and fn+(^~^ = f'^'~°'~^(,Un+a, ■ ■ -.Un+Ai-p) is a constant 
function or such that 

r.rM-a-13 r.rM-a-13 
J n+a ^ J n+a _^ g 

dUn+M-p 

for some integers a^j3 > 1. 

Performing a local change of the variables Un i — > Un = f{un), we can reduce to 

gKliUni • • • 5 Un+Kl) = f^^ ° ^{Un+a^ • • • 7 ^n+A/-/?)- 

Finally, we prove a = p. Indeed, on the one hand, from equation 

[M, M - /3] {{Un, Un+Al} , Un+2M-f}} = {«n, {Un+M, Un+2M-I3}} 
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after some elementary computations, we have 

^Og gpiiUn+a, ■ ■ ■ , Un+M-l3),u„+M-p 
II 

gM-l3{Un+M, ■ ■ ■ , Un+2M-(5),u„+MgMiUn+M-l3+a, ■ ■ ■ , Un+2{M-I3))~^ 

where the variables appearing on the left-hand side do not intersect with those ap- 
pearing on the right-hand side. Then, there exists a non-zero constant k, such that 

QM-piUn+M-, ■ ■ ■ , Un+2M-I3),u„+M = k gM{Un+M-l3+a, ■ ■ ■ , Un+2iM-l3)) 

and this PDE makes sense only if a > /3. 

On the other hand, repeating the same argumentations for equation 

[M - a, M] {{Un, Un+M-a} , Un+2M-a} = {Un, {Un+M-a, Un+2M-a}} 

that is 

gM-a{Un, ■ ■ ■ 1 Un+M-a) ,Un+M-a9M{Un+a, ■ ■ ■ i Un+M-^) 

II 

log A/ (Wn+M, • • • , Un+2M-l3-a),u„+M 

we arrive aX a < (3, that implies a = (3. □ 

For any fixed order M, the leading order functions, given by Lemma 2.2, define 
essentially different classes of PBs (2.1). 

Definition 2.3. Let (a,^) be a pair of non-negative integers with a > 1, we call 
(a, i^)-brackets the class of non-constant PBs (2.1) of order M = 2a + ^ expressed in 
the canonical coordinates, i.e. guiun, • • • , Un+ivi) = f^iun+a, ■ ■ ■ , Un+a+^)- 

For the constant case / = cr^/ one can immediately prove the following 

Proposition 2.4. // there exists a set of coordinates reducing the leading order term 
to the constant form gM{un, ■ ■ ■ , Un+Ki) = <Jm, for some non-zero constant aM, then 
all the coefficients gk, k = 1, . . . , M are constant in such coordinates, i.e. 

{Un, Un+k}M = 9k{Un, • • • , ^n+fc) = CTfc, k = 1,...,M, 

where ak are complex constants. 
2.2. Classification theorem. 

Theorem 2.5. For any PB (2.1) of order M , there exist a set of coordinates and an 
integer a > 1, such that the coefficients gk{un, • • • , Un+k), k = 1, . . . , M are given by 
linear combination of the suitably shifted function f^ (see Lemma 2.2) 




where {Ti}-^Q ^ are constrained complex parameters. Explicitly, we have 

(min(a+p,a) \ 
s=max(0,p) / 

where ^ = M — 2a, p = —a, . . . ,a, T is the shift operator (1.3) and Af„^ are scalars 
such that 

(i) they satisfy the multiplication rule 
(2.7) K-r = K-r K-s s,r = 0,...,a, 

(ii) denoting 9 = min(Q!,^), they can be expressed by explicit formulas in terms of 
the 9 + 1 parameters Tq, Ti, . . . , Tg^i, 



(2.8) 



K 

A!^ 



det Ar ({(-) Vq-^ t,}s>o) r = 0, . . . , a - 1 



^0 (.n) 



where Ar is the band- Toeplitz matrix 

^ a.\ ao 



(2.9) 



Ar (a) 



\ 



V 



«r-l 
Or 



; ; ; ■■. 

... (1\ CIq 

. ... 02 ai J 

associated to the sequence a = (oq, ai, . . .). 
Here and below we adopt the notation Tp = if p < or p > 

Remark 2.6. The scalars appearing on formula (2.6) can be easily visualized looking 
at the rows of the following rhombus 



\a—l \a 











... A^^ 


A? 


(2.10) 




K 


AS 


\ a—l 

... Aq 


A? 






A", 


Ail 


... A^I^ 


A^l' 



A° 
A". 



Q + l 



A^ 



-a+l 



We complete the description above, specifying the non-trivial constraints that the 
Jacobi identity imposes on parameters r in the following 

Theorem 2.7. Let (a,^) be a pair of non-negative integers and define the sequence 
{(Jq}q>o = {(— )'^A"_g}g>o. Then, the [a, brackets are in one-to-one correspondence 



(2.13) (ror+Met^5-p({^J</>o) = (r^^+Met ({a„_ J, 



with the intersection points of the 9 projective hypersurfaces in CP defined by the 
homogeneous polynomial equations 

(i) if ^ ^ — 1 , for any p = 0, . . . ,a — 1, 

(2.11) (ro)P+Met^5_p({aJ,>o) = {T^y+' detApi{a^^,},>o) 
(a) if ^ < a, for any p = 0, . . . - 1, 

(2.12) (-ro)^"% = (T5)«-^+Met^5„p({a„_J,>o) 
(Hi) if a < ^, < 2a — 1 , for any p = 0, . . . ,^ — a 

}g>Oj 

and, for any p = ^ — a + l,...,Q; — 1 
(2.14) [r^T^r^ = in)^""^' det A^.p ({a«- J,>o) • 

The case 6' = is described below in the Example 2.4.1. 

2.3. Proof. In our computations we often encounter generalized Fibonacci sequences. 
Therefore, it might be useful to remind some elementary topics about them. 

Remark 2.8. For any positive integer k > 2, the /^-generalized Fibonacci sequence 
{Fn}n>o is defined by the order k linear homogeneous recurrence relation 

K = Eli(-r'«^i'n-., forr2>l, 

for arbitrary coefficients ai, i = 0, . . . , k. 

The generating function for {F„}„>o is given by 

n>0 

and, using the theory of lower triangular Toeplitz matrices (see for example [13] or 
[16]), one can express the /c-generalized Fibonacci numbers in terms of determinants 
of matrices with entries given by the coefficients of their recurrence equation, that is 

(2.16) F„ = det(ai_i+j)i<ij<„ = det (a) , 

where a = (oq = l,ai, . . . ,ak, ak+i = 0, . . .) and A the band Toeplitz matrix defined 
by (2.9). 

Let us consider a PB of type (2.1). According to the Lemma 2.2 we can choose certain 
coordinates on the target manifold Ai'^ such that the leading term gM reduces to the 
following form f^+a = f^{un+ai ■ ■ ■ :Un+a+()- OUT opiuion, the simplest way to 
obtain the form of coefficients g^, k = 1, . . . , M comes from considering equations 

[M — a—p, M] {{Un, Un+M-a-p} , Un+2M-a-p] = {^ni {Un+M-a-p: Un+2M-a-p}} 
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[M,M - a - p] {{Un, Uri+M] , Un+2M~a~p} = {Un, {Un+M, Un+2M~a~p}} 

for any p = 0, . . . , a + ^ — 1. 

2.3.1. Formula (2.6). We start a detailed analysis of some initial cases (i.e. p = 
0, 1), that will suggest us how to organize the general computations. For the sake of 
simplicity, we suppose ^ > a. 

p = 0. Equations [M - a, M] and [M, M-a] yield (log/|) = tq and (log/«) 
for some non-zero constants tq, (see condition (2.4)) and 

where dots stay for any arbitrary function depending at the most on the variables 

Un+i, and A[j = (to)"^ r^, Ag = Tq (t^)"^ 

p = 1. Analogously, equations [M — a — 1,M] and [M,M — a — 1] provide us 
{^osfn)^u„+i " ^1 = T-^-i for some constants Ti, r^^i and 

5'a+5— l('Wn) • • • ) l),M„+Q+5-l 
II 

T'O fi'Af-l('Un, • • • , Un+2a+^-l) + ^1 f^{Un+a, • • • , 'Wn+a+5) 

5'a+5-l('f^n, • • • , Mn+a+§-l),-u„ 
II 

It necessarily follows that the constants Ti, r^_i are non-zero and 
with X'^zl = iny^n-i and A-_i = -(ro)"iri. 

Remark 2.9. Notice that the function hf^^ can be understood as the leading order 
function of a (a, ^ — l)-bracket: /i^+l^, = f^~^{un+a, ■ ■ ■ , Un+M-a-i)- 
We decide to forget about the contributions provided by the leading order functions 
f^' of lower order PBs (i.e. M' < M), postponing to Section 4 the problem of 
classifying the compatible pairs of PBs (2.1). 

According to Remark 2.9, we solve equations [M — a — 1, M] and [M, M — a — 1], 
finding 

ga+^-l{Un, ■ ■ ■ ,M„+a+g-l) = \\ /| + . . . + A^^ Vn+a-l' 

where A°^ = (Aq)~^A"„i, A°7^ = ^^o -^ali and dots stay for any arbitrary function 
depending at the most on variables Un+i, • • • , Un+a+^-2- Iterating this procedure, one 
can show that the coefficients {gk}k=i,...,M depend on lattice variables according to 
formulas 

9a+^+k = ga+S,+k{Un+k, • • • ; Un+a+() i k = 0, . . . , a 

Qa+i-k = ga+^~k{Un, . . . , Un+a+^-k) , k = 1 , . . . , a + ^ — 1 
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and 



where 



, Un+^) = A° exp Ti Un+i^ 



(i) A" is a free multiplicative constant, that can be normalized (i.e. = 1) 
choosing a suitable rescaling of the coordinates: u„ i — y k Un, for some constant 
k, 

(ii) {Tj}-^Q ^ are non-zero complex parameters, thanks to condition (2.4). 
We denote Tj = 0, if z < or 2 > ^. 

Moreover, solving the equations [M — a — p, M] and the symmetric ones [M, M — a — p], 
we obtain the formulas (2.6) for coefficients {gk}k=i....,M- 

2.3.2. The multiplication rule for constants X's. Let us focus our attention on equa- 
tions [M — a — p, Af]p^Q a+^-i (analogous results come from equations [M, M — a — p]). 
The constraints on constants A's can be encoded into the following linear system 



/ r, XT/ \ 




(2.17) 



v 







I 



A 



a—p+l 
a—p 



\ 



\ Q— 1 

^a—p 





> a—p+l 



\ a— 1 

^cv-p+1 

\a 

^a-p+1 



a-1 
a-1 



a-1 



\ 





K ) 



( ro \ 

n 

Tp~l 

\ rp / 



where Af_^ is non-zero iff s, r = 0, . . . , a, and recursively on p we have defined 



(i) Xlp^ — Aq A"_p, 

(ii) A^lf^ : 



Combining (i) and (ii), by induction on p one can prove the following multiplication 
rule 



(2.18) 



A! 



S, r = 0, 



, a, 



that allows us to describe all constants A's in terms of the constants {A"_^, XTVis=o,...,a, 
appearing on the edges of rhombus (2.10). 

2.3.3. Constants X 's as function of parameters r 's. We start looking at the last row of 
the matrix appearing on (2.17). If < p < a, we have the following order p recursive 
relations 

(2.19) XTp=-iro)-'i2^^XTp+^■ 



Denoting Fp = A"_p, Up = (— )^t"q ^ Tp, we recognize the p-generalized Fibonacci se- 
quence (see equation (2.15)). Then, according to Remark 2.8, we can express 

K-p {'To, ■■■,rp) = det Ap {{{-Ytq^ r,},>o) p = 0, . . . , a - 1. 
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To write down the constants {XaZl}s=o,...,a as functions of the constants {\a-s}s=o,...,a, 
we look at the equations 

[M - q,M - a + g]g=0,...,a {{Un, Un+M-q} , Un+2M-a} = {Un, {Un+M-q, Un+2M-a}} ■ 

They spht into 

(b) g(a+0+q i'^n+q, ■ ■ ■ , Un+a+i)) ,Un+q = fK'^n+q, ■ ■ ■ , Un+q+^)) ,u„+q 

that are compatible only if 

(2.20) A^ = (A^)"^A^,. 

Substituting (2.20) in the formula (2.18), we hold the multiplication rule (2.7). 

2.3.4. Constraints for parameters t's. Supposing ^ > a, we can analyze when a < 
p < ^. Analogously to (2.19), we have the following linear homogeneous recurrence 
relation of order a 

a 

(2.21) Tp = - ^ K-qTp~q^ fo^ ^ny p>a. 

9=1 

For any q = 1, . . . ,a, denoting a, = (-)'^A"„g, we arrive at Tp = Yl'l=ii.~Y^^(^q^P-q 
that is 

(2.22) Tp(ro,...,r«_i;r^) = ro det^p({(-)'A:^_J,>o) , p = a,...,C. 

Finally, when ^ + 1 <jo<a; + f, we obtain analogue relations with respect to the 
case 1 <p < a: the equations [M — a — M] give us 

q 

= 0' g = !,...,« -1. 

Following again Remark 2.8, for any g = l,...,a — 1, we obtain 

(2.23) T^.q (to, . . . , r„_i; r^) = det Aq ({(-)^(A^)-i A:},>o) 

Observe that formulas (2.22) and (2.23) have to be consistent: this provides us the 
constraints on parameters r's. According to the values of a and ^, we immediately 
find out the homogeneous polynomial equations given by Theorem 2.7. 

2.3.5. Remaining equations. Due to the relation (2.20), the symmetric set of equa- 
tions [M, M — a — p]p=o a+^-i gives us the same contraints. 

Moreover, replacing expression (2.6) on the other equations [p, gj^ j^^, with straight- 
forward computations that generalize the previous ones, we can obtain that the rela- 
tions for the parameters tq, . . . , Te_i; describe above are necessary and sufficient. 

2.4. Examples. Looking at the number of the lattice variables which the leading 
order function can depend on, we study in detail the minimal and the maximal cases. 
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2.4.1. (a, 0)-brackets. The leading order coefficient can be choosen of the form = 
f^{un+a) = exp(roti„+Q) and substituting in the formulas for the constants A's we 
obtain 

{Un, Un+a} m = ^Xp (tq M„) + exp (tq Un+a) 

where Tq is a free complex constant. Notice that the cubic PB of Volterra lattice (2.3) 
belongs on the class of (1, 0)-brackets. 

2.4.2. {1, brackets. These PBs are given in the canonical coordinates by the for- 
mulas 

{Un,Un+M} = exp(2;„+i) 
{Un,Un+M-l} = TQ^T^exp{Zn)+ToT^'^exp{Zn+l) 
{Un,Un+M-2} = exp{Zn) 

where Zn = Yli=o(~^y^o Un+i and the pair (tq, r^) belongs to the set of points of CP^, 
described by the following equation (see equation (2.11)) 

(2.24) ro"+^ + (_)"+€r^-+€ = O where [tq : t^] G CP\ « = 1. 

Analogously, the (a, l)-brackets are parametrized by the same points (2.24), with 

2.4.3. A multi-parameters example: {2 .,2) -brackets. This is the ffist non trivial case 
in which the parameters r's are described by the intersection points of certain hy- 
persufaces. From the Theorems 2.5 and 2.7, it follows that the (2, 2)-brackets are 
characterized by three parameters To,ri,r2, constrained to satisfy the following sys- 
tems 

{r,f = {r^f 

3. A Darboux-type theorem 

In Theorem 2.5 we have proven that, up to local point-wise change of coordinates, 
the leading order function qm of any PB (1.5) can be considered of the following form 
9m = fi+a = exp (Yli^Q Ti Un+a+i^ , for somc pair of non-negative integers (a, ^). 
In the present Section, we deal with the problem of reduction of PBs (2.1) to a 
canonical form by more general changes of variable than the local ones. Let us give 
the following 

Definition 3.1. Let (si,S2) be a pair of non negative integers, we call a discrete 
Miura-type transformation any map of the form z„ = Lp{un-si, ■ ■ ■ i'^n+s2) that is a 
canonical transformation, i.e. a change of coordinates preserving the PBs. 

Notice that these discrete Miura-type transformations, as the continuous ones, are 
differential substitutions (i.e. they depend on Un+i,Un+2, ■ ■ ■) and therefore, only 
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formally invertible. Defined new coordinates according to 

(3.1) Zn = Tj Un+i , 

we prove by direct computation the following 

Theorem 3.2. Any {a, ^)-bracket is mapped by the Miura-type transformation (3.1) 
into the following [a + S,-, 0)-bracket, 

/g 2) {Zn, Zn+2ia+o} = Tq exp(z„,+a+^) 

{Zn,Zn+a+^} = Tq [exp{Zn) + exp{Zn+a+i)] ■ 

Remark 3.3. Subdividing all the particles- variables into a + ^ families, according to 

= 2:{„+^){n-i)+p, p = l,...,a + ^, 
the PB (3.2) can be presented in the following simple form 

v'n , r = ^on exp [V'^l^ ) 



(p) (p) . 



Therefore, it splits into a + ^ copies of the cubic Volterra PB (2.3). 

Furthermore, any constant PB of order M, Wn+fc} = (^ki with A; = 1, . . . , M, can 

be reduced by the lattice splitting Vn'^ = UM{n-i)+p, to the following PB 



(3.3) 



,,(p) ,,(9) . - 

(P) (<?) \ _ 



where ctq = and p, 5 = 1, . . . , M. Notice that, if all constants are normalized (i.e. 
CTfc = 1, A; = 1, . . . , M), the PB (3.3) becomes the quadratic PB for the Bogoyavlensky 
lattice (BL) of order M (see the definition of BL in Suris [12], Chapter 17). 

Proof. We are interested on the expression of (a, ^)-brackets in z-coordinates. Let be 
p = 0, 1, . . ., we have to compute 

= J2i=0 '^i J2j(^J+(i,p) T?' {""n+i; ^n+p+j} + J2j£,J-{i,p) '^j {^n+p+j, Un+i} 

where we have defined the sets of admissible j's, 

J+{i,p) = {0,...On{z + ^-p,...,z + 2a + e-p} 
J-(z, p) = {0, . . . O n {z - 2a - e - P, . . . , ^ - e - p}. 

In the following, to avoid some technicalities, we suppose ^ > a. At first, we notice 
that when p > 2(a + ^), {zn, Zn+p} vanishes. 

Let us start considering the set J~{i,p). It is non-empty only if i = ,^ and p = 0. 



15 



When p = 0, J^{i,p) is non-empty only if z = 0, then our summation vanishes, indeed 
If p > 1 , we can reduce to evaluate 

(3.4) {Zn, Zn+p} = X]i=0 '^jeJ+{i,p)^j {'^n+i,Un+p+j} ■ 

In the following, we give some details about the complete computation. Enforcing 
a recursive procedure on p, for any i that runs from to ^, we describe the set of 
admissible j: J^{i,p). 

Step 1: p = 2{a + ^). Then J^{i,p) 7^ if and only if z = ^ and 

Step 2: p = 2{a + ^) — q, q = 1, . . . In the following table: fixed p, we describe 
the non-empty sets J^{i,p), as the index i changes. 



p i: J+ 


{t,p) ^ 


Jtihp) 






2(a + 0-g i = i 




{q-2a 


• g-1 


?} 






{ 




...} 




-q + 2a 


{ 1 ... . 


. 2a -1 


2a} 






{ ■-. ... . 




...} 


1 = ^ 


- q + a 


{ . 


a — 1 


a] 






{ 




...} 


1 = ^ 


-q+1 


{ 





1} 


i = i 


- Q 


{ 




0} 



Looking at sets J^{i,p), we distinguish three cases depending on the number of 
elements ^J^{i,p): 

(^) i^J^hp) = 2a + l 

(n) 1 < #J+(2,p) < a + 1 
(in) a + 1 < ij^J^{i,p) < 2a + l 

(i) When j^J^{i,p) = 2a + 1, we have 

^ ] ^ ] Tj {Un+i^Un+p+j} = '^^Ti '^^Tr-t{Un+i,Un+i+M-t} 

i=0 je{i+^~p,...,i+2a+^-p} i=0 t=0 

where r = i + M — p. Now, 

s=0 -^a-sJ I'^n+a-S) • . . ; Un+a+^-s) Z^t=s '^''-t^a-t+s 

and Yln='^ '^T-tKx-t+s — Y11=Q'^r~s~qK>i-q = 0' according to the recurrence relations 
(2.21) for the parameters r's. 
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(ii) When 1 < i^JaihP) < + 1, the set J^{i,p) is given by J^{i,p) = {0, . . . ,r} 
for some < r < a. The summation (3.4) can be written in the following way 



i=^-q 



je{o,...,r} 



r = i — ^ + q. 



Noticing that Yl '^j {Un+UUn+p+j} = Y rr-s{Un+i,Un+i+M-s} = To Xa-rfn+i+a- 
j6{0,...,r-} s=0 

we obtain 



^-q+a 
z=^-q 



j€{0,...,r} 



'^-q+a 



J 7 



and Yi=l~q '^i K~^+q = Yl'i=o^-q+tK = 0' ^cc the rccurrcnce (2.21). 
(iii) When a + 1 < jj^J'^ii.p) < 2a + 1, the set Jaihp) is JitihP) = {0,...,r}, 
for some a < r < 2a. With analogous calculations, one can directly prove that the 
summation Yje{o,...,r}^j {'^n+i,Un+p+j} vanishes. 

Step 3: p = 2a + ^ — q, q = 1, . . . ,a + 1. The sets J^{i,p) are given by 



P i-Jti.hP)i^0 






2a + i-q i = i 


{^-2a + q . 


■ e } 




{ ... 


. ... ••. } 




{ ... 


■•• } 


i = i-q 


{ e-2a . 


e-i 




{ ... 


} 


i = 2a — q 


{ 


2a} 




{ 


■ } 




{ 


■■■ } 


% = 


{ 


... q} 



We find out again the two situations 

(z) i^J^i^,p) = 2a + l / 

(a) a+l < i^J+{i,p) < 2a + 1 / 

that we have been already studied in the previous Step. Then the summation (3.4) 
vanishes. 
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Step 4: p = a + ^. In this case, we have 













a + ^ {i — a, . 


+ a} i = i 




. ... ^ 


} 






{ ... . 




} 






{ ... . 




■■• } 




i = ^-a 


{e-2a . 




... 










} 




i = a 


{ . 




2a} 






{ 




} 






{ 




} 




i = 


{ 


... 


a} 



three different cases 

(i) = 2a + l / 

(a) a + l < ^J+{i,p) < 2a + 1 / 
(m) = « + l ? 

Only the last one provides some contributions. Indeed, 

and 

^j<:{(-a,...,(}'''j {'^n+i,Un+p+j} = Xq [Un+a+i, ■ ■ ■ , Un+i+a+£,) ■ 

Our summation (3.4) becomes 

{Zn, Zn+a^} = ^0 ^ f{Zn) + f (Zn+a+i) = Tq [exp{Zn) + exp^Zn+a^)] ■ 

Finally, Step 5: p = 0, + a — 1 can be analyzed similarly to Step 2 and Step 

3. □ 



4. Compatible pairs 

We are now in a position to study the compatible pairs (P, P') of PBs (2.1). We pro- 
vide some necessary conditions that we expect to be also sufficient. This might be the 
starting point for a future classification of the still little-understood bi-hamiltonian 
higher order scalar-valued difference equations. 

Let P and P' two PBs of type (2.1). Their leading order functions are given respec- 
tively by 

(4.1) 9M{Un, Uu+m) = akliUn) fH^n+a, ^n+a+c) (Wn+A/) 



(4-2) g'M'i^n, • • • , Un+M') = a\p{Un) /^'{u^+a', ttn+a'+f ) ^^/'(^n+A/') 
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where /^+^ = aj,/ exp (Ep=o and /^^^, = ciA/'exp (Ep=o ^n+a'+pj 

for some non-zero constant aM and a ml It is not restrictive to suppose M > M' . 

Lemma 4.1. A pair of non- constant PBs (2.1) {P, P') forms a pencil of PBs only 
if there exists a local change of variables, reducing the leading coefficients (4.1) and 

(4.2) to the formulas 

(4.3) giVliUn: . . . , Uu+m) = f^{Un+a, • • • , Un+a+^) 

and 

(4.4) g'^j, {Un,..., Un+M' ) = (Wn+a , • • • , ) 

for some functions f^ and f^'. Notice that a = a' . 

Proof. Let P and P' PBs defined respectively by leading function (4.1) and (4.2). 
From equation 

[M, M'] {{Un, Un+M'} , Un+M+M'} = {Un+M', Un+M+M'}} 

we find out (log aM(^in)),„„ = (log a'^j, («„,)) , that implies aM{un) = ka'^piun) for 
some constant k. A suitable change of coordinates leads us to (4.3) and (4.4), where 
the constant k has been absorbed in the multiplicative constant (Tm'- We hold a = a' 
looking, for example, at equation 

[M, M' - a] {{Un-, Uu+m} , Un+M+M'^a} = {Un, {Un+M, Un+M+M'-a}} 

that makes sense if and only if M' — a < M' — a\ i.e. a > a'. In such case we have 

(4-5) •^"+",M„+Q+5 fri+a+£,+a' = 9a+('i'^n+M, ■ ■ ■ i Un+a+^+M') ,u„+m fn+a 

that gives g'ci+^i{Un+M^ ■ ■ ■ , Un+M+M'-a) ,Un+M = ^ (Un+a+^+a' , • • • i Un+a+a' +^+^') ■ 

This equation makes sense if a < a'. It follows a = a', so it is not restrictive to 
suppose C > n 

Theorem 4.2. A pair of non- constant PBs (2.1) (P, P'), defined by leading functions 
(4.3) and (4.4) for certain parameters r and r' satisfying algebraic constraints of 
Theorem 2. 1, forms a pencil of PBs only if 

(4-6) = 7p = T-^-^'+p, 

for any p = 0, . . . 

Proof. Let P and P' be a pair of PBs (2.1) respectively, of order M = 2a + ^ and 
M' = 2a + ^' (i.e. M > M'), leading functions (4.3) and (4.4) and suppose that any 
linear combination P + u P', for /i, u arbitrary constants, is a PB of order M (i.e. 
it satisfies the bi-linear PDEs, coming from Jacobi identity). 

Analogously to the procedure followed in the proof of Theorem 2.5, we find necessary 
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conditions looking recursively at equations 

[M', M -a-p], [M - a - p, M'] p = 0,...,a + ^' 
When p = 0, equation 

[M', M — a] {{Un, Un+M'} , Un+M+M'-a} = {Un, {Un+M' , Un+M+M' -a}} 

gives us 

log/^ {Un+a, • • • , Un+aH') ,u„+^+i., = K ^'^SfH'^n+M', • • • , 'Un+A^'+^ ) 

that enable us to identify r^, = r^. Analogously, equation[M — a, M'] provides Tq = Tq. 
Iterating this procedure, when 1 < p < a, equations [M — a — p, M'] restricted to the 
function originate some constraints that can be organized in the following matrix 
form 







V 



/ x/a-p 
^ct—p 



J 



p+1 \/Q— p+1 
'^a—p ^a— p+1 



v 



p 

^a—p 



\ to— 1 
^a-p+1 
\ to 

^a-p+1 



\ to-1 





A'" 



a / 



Tp-l 



and recursively on p, we obtain that A^_p = A^.^, for any p = 0, . . . , a. Adding the 
contribution coming from equations [M', M — a — p], the following constraints on the 
parameters r's hold 



n-p - '^/-p 

The formula (4.6) immediately follows. 



p = 0, . . . , a — 1. 



□ 



4.1. Compatibility with constant brackets. We devote this sub-section to study 
the compatibility of pairs (P, P') of PBs (2.1), where P' can be reduced by local 
change of coordinates to the constant form. 

Lemma 4.3. Any non-constant PB (2.1) of order M = 2a + ^, with C. > is 
compatible with a constant bracket of order M' only if M' < a. Moreover, the constant 
coefficients {o"j}j=i,,,,,A/' have to satisfy the following recurrences 



(4.7) 



0. 



Proof. Let us fix a pair of PBs (2.1) (P, P') such that their leading functions are of 
the form 

dAliUn, ■ ■ ■ , Un+Nl) = 0,M{Un) f^{Un+a, ■ ■ ■ , ^n+a+c) OA^(^^n+A/) 
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As in the proof of Lemma 4.3, from equation [M, M'], we obtain (log aj\f ^ = 
(log a^j,(w„)) that is aAf(wn) = kaM'{un), for some constant k. 
Then, after a change of variables, we can reduce to consider 

guiUn, ■ ■ ■ , Un+u) = f^{Un+a, • • • , Un+a+£,) 
g'wi'^n-, • • • , Un+M') = M' ■ 

Let us now suppose a < M' < M, then from the equation 

[M, M' - a] {{Un, M„+m} , M„+M+M'-a} = {Un, {Un+M, Un+M+M'-a}} 

we obtain f^{un+a, • • • , ^in+«+«),«„+„+^o-M' = that implies cta/' = 0. 

Let us suppose > and k = 1, . . . , a, then M — k > a. We focus our attention on 

equations 

[k, M - k] {{Un, Un+k} , Un+Al} = {Un, {Un+k, Un+Al}} + Uu+m} , Mn+fc} , 

obtaining the constraints system A^I^cTq, = Tq^ Yli=o ^a+i-k that can be written in 
the form (4.7). 

□ 

Remark 4.4. When ^ = 0, we immediately find that any (a, 0) -bracket is compatible 
with a constant PB of order a, given by the formulas 

{Un,Un+a-s} = 0, S = 1, 1. 

4.2. Examples. We complete this Section, adding the details for the relevant family 
of (1, -brackets, where we are able to prove that our necessary conditions are also 
sufficient. Let P be a (1, ^)-bracket, 

- if To 7^ r^, according to Theorem 4.2, there are not (1, -brackets P', such 
that the pair [P, P') forms a pencil of PBs. Looking for constant brackets, we 
have that P is compatible only with first order constant brackets if and only 
ife = 0,L 

- if To = Tj, then necessarily M = 2K for some positive integer K. Any pair of 
(1, 2q')-brackets, with q = 1, . . . , K, defined by building functions of the form 

where Zn"^^ = Yl'i=o'^\~'^o)'"^n+i7 forms a pencil. 

5. On non-degenerate, vector-valued PBs 

In order to provide some examples of vector-valued PBs, we can define new lattice 
variables, according to the following 
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Proposition 5.1. Let us consider a pair of positive integers {M,K), with K < M. 
Any scalar-valued PB (2.1) of order M, according to the formula 

^n"^^ = ^ni^+p> p = 0,...,K -1, 

is transformed into a non-degenerate PB (i.e. the leading order is given by a non- 
singular matrix) of order A and target space of dim. K iff M = AK, for some 
positive integer A. 

We are interested on vector- valued PBs of first order (i.e. A = 1), tlien we put K = M 
in tlie above proposition. First, let us recall some preliminary definitions. 
A Lie group G, with a Poisson bracket { ■ , ■ j^, is a Lie-Poisson group if the mul- 
tiplication yUiCxG— T-Gisa mapping of Poisson manifolds, where on G x G is 
defined the bracket {ip, V'IgxG i9,h) = {ip{,h) ,1^ {, /i)}^ (g) + {^p {g , ) , ip {g, )}^ (/i), 
with h,g & G. Let cf^ be the structure constants of a Lie algebra g. The couple (g, 7) 
is a Lie bi-algebra if and only if 

(i) 7 is a 1-cocycle on g with values on g ® g, where g acts on g ® g by the adjoint 

(2) 

representation ad^ = adg ® 1 + 1 (8> ad^, that means: (5ad7 = 0, i.e. 
adf (7(^))-adf (7(O)-7([e,^]) = 0, 

or, fixed a basis of g, €^^7^? = €^^7^'' + c^s7r ~ '^^rll'^ ~ (^IrlT- 

(ii) *7 • S* ® S* — ^ 0* defines a Lie bracket on g*: 77]^. =* 7 (.^ r]) . 

The correspondence between the Lie-Poisson groups and the Lie bi-algebras is clarified 
by the following 

Theorem 5.2. Let G a Lie group, with tangent Lie algebra q: locally a Lie-Poisson 
structure on G is uniquely (up to isomorphism) determined by a Lie algebra structure 
on the dual space g*, then q is a Lie bialgebra (g,7). 

Proof. All the detail of the proof can be found in [8] . Here we recall some ideas about 
the direction from Lie bi-algebra to Lie-Poisson group. Fixed a basis of (g,7), the 
constants 7^*^ define a Lie bracket on g*. Moreover, solving the differential equation 

_ f),TrP1\ 

Ik — '^kT^G le' 

we can define = dp(j)'K'Q dqip , that is a PBs because the compatibility condition 

of the system 

dkT^G = C^fe^G + + 7r 

is guaranteed by [i), i.e. 7 is a 1— cocycle on g. □ 

We focus on the sub-class of Lie Poison group, given by the following 

Definition 5.3. A Lie-Poisson group {G, { , }c} is called admissible, if there exist: 
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(i) A skew-symmetric matrix k6 A^g, such that the cohomologous 1-cocycle 7, 
defined by 7 = (^adk + 7, i.e. 7^ = 7^^ + c^^/c"*^ + k^^^cl^, provides a Lie algebra 
structure on q*. Notice that k has to satisfy the Yang-Baxter equation 

(ii) A Lie algebra homomorphism r : (0*,7s^) — > idiCpg) such that r* : rj'^ = r^^ 
defines a Lie algebra homomorphism : (0*,7f^) (0,0^^). 

We are now in a position to formulate the Dubrovin's theorem [3]. 

Theorem 5.4. An admissible Lie-Poisson group [G, { , }q) together with correspond- 
ing matrices r, k defines a Poisson bracket of the form 

.52) {<><}i = ^'Mun) 

{<,<+i}i = ^'-^ (u„,u„+i) 

where u„ G G for all n, according to the following formulas 

{if (Un) , (u„+i)}^ = da^P (U„) r"^(9^?/) (u„+i) 

where da and d'^ are left- and right-invariant vector fields on G, 

{(f (u„) , ifj (u„)}^ = h"^ (u„) daip (u„) dfji/j (u„) 

where /i"'^(u„) = tTq^ (u„) + Ad|^^2i/c°^ and 7r^^(u„) zs determined by the system (5.1). 
i/ere y?, "0 are arbitrary smooth functions on G. Viceversa, all brackets of such form 
are obtained in this way under the non- degeneracy condition det g^^ 7^ 0. 

Remark 5.5. This theorem does not seem to have simple applications in the lattice 
systems, studied in the literature. For example, one can notice that the fundamental 
Toda lattice has three well-known local compatible PBs (see [12]), but only the qua- 
dratic one is non-degenerate and can be represent using the previous theorem on the 
algebra of Aff°M^, group of affine transformations of the straight line. 

Our new vector-valued (A^ > 1) PBs (5.2) are provided by the following 

Theorem 5.6. Any {a, S,) -bracket, after the consolidation lattice procedure, becomes 
a non- degenerate PB of the form (5.2), with associated Lie bi-algebra given by 

0(a,O = span {^s},=i M ' [Lp. ^9] (^n) = C^g Lr(v„,) 

S*,^^) = span{/2^},^, ...^,,, [R^,R^] (v^ = 7^ ^^(vn), 

where the summation is over the repeated index r and c^^ and 7^'^ are certain constants 
depending on the parameters r 's. 

Proof. Choosing A' = M in the Proposition 5.1, we have 
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and h^^{\n) = 9'^^{^n,^n) ~ fi'*''(v„,; v„). Substituting the explicit formulas (2.6), we 
find out a decomposition of the leading order (7*'' (v„, v„+i) = L^(v„)i?^-' (v„+i). The 
matrix L(v„) is given by L(v„,) = diag(/i, . . . , /Af)(v„) where • • • , <"^^^') = 

exp (eLo '^i and 



Aj 



AS 



V AO ... X-J 

Analogously the matrix R{-Vn+i) f actor izes into i?(v„_|_i) = diag(ri, . . . , r^) A^, where 



") = exp(Er=o-^v,_,<;r), 

/ A" 



A 



R 



\ 



\ ... A^/ 

Notice that we are using the notation = if s < 1 or s > M. Moreover, to avoid 
some technicalities, we suppose ^ > a. 

Now, let {Ls}s=i,...,M be the s-th column of matrix L and Q{a,^) = span {Lg} as vector 
space. By direct computation, we find that the following commutators 



(5.3) 



[Lp, Lq 



equip 0(o,5) of a Lie algebra structure. At first, we observe that if a + < < M, 



Lq{wn) are constant vector fields. Then 



p=l,...,M, 



for some constants apg that can be expressed in terms of parameters r. 

In particular, observe that [L^, La+^+t]'^ (v„) = 6s.t Ls, where s,t = 1, . . . , a and 6 is 
the Kronecker symbol. 

Let us now consider p,q = 1, • • • , Q.+^. It is not restrictive to suppose q > p. Denoting 
t = q — p, the formulas (5.3) for the commutators reduce to 



[Lp, L. 



p+tl 



We distinguish two cases: 
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(i) When t < a, we have [Lp, Lp+tf (vn) = Xa~ktl XlLo W(vn)^p(v„). 
Let us focus on the summation Yll=o K-l'^i- According to formulas (2.7) and 
(2.21), we have 

t t 

i=0 i=0 

(ii) When t > a, we obtain [Lp, Lp+tf (v„) = X'^ZlXl ^^"=0 K-l^t-a+i lp+t{yn)lp{yn) 
and we have Yl'^^oK'^t-a+i = K Y.'^=oK'^t-a+i = 0, according to (2.21). 

Finally, starting from the rows of matrix R, we can repeat the procedure above finding 
the Lie algebra structure on g*. □ 
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